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Abstract 

We consider hadron production in deuteron gold collisions at RHIC in the forward 
rapidity region. Treating the target nucleus as a Color Glass Condensate and the 
projectile deuteron as a dilute system of valence quarks, we obtain good agreement 
with the BRAHMS minimum bias data on charged hadron production in the for- 
ward rapidity (y = 3.2) and low pt region. We provide predictions for neutral pion 
production in minimum bias deuteron gold collisions in the forward rapidity region, 
y = 3.8, measured by the STAR collaboration at RHIC. 



1 Introduction 



The recent observation of the suppression of the charged hadron spectra in the forward 
rapidity region P at RHIC, and its centrahty dependence has generated a lot of interest and 
excitement in the high energy heavy ion community. The suppression of particle spectra 
and disappearance of the Cronin effect at forward rapidity had been predicted ^ |3J in the 
Color Glass Condensate formalism |;4j, unlike the more conventional models which predicted 
a stronger enhancement of the spectra in the forward rapidity as compared with mid 
rapidity jH] • The forward rapidity data strongly suggest that the high gluon density region 
of QCD phase space, the Color Glass Condensate, has been observed at RHIC. While the 
Color Glass Condensate has been successful in predicting some global features of the data, 
such as multiplicities and their energy, rapidity and centrality dependence in gold gold and 
deuteron gold collisions 6J as well as qualitative predictions [21 E] of the suppression of the 
hadron spectra and their centrality dependence, there has not been a quantitative analysis 
of the forward rapidity hadron spectra using the Color Glass Condensate formalism. Here, 
for the first time, we provide a limited, but quantitative analysis of the low pt, forward 
rapidity RHIC data, using the Color Glass Condensate formalism. 

As emphasized in j7j, the forward rapidity region at RHIC is the best kinematic region 
to look for the signatures of the Color Glass Condensate since this is the region where one 
probes the smallest x in the target nucleus so that the Color Glass Condensate will be 
manifest more strongly in this kinematic region. Also, forward rapidity deuteron gold col- 
lisions are ideal since one does not have the final state (Quark Gluon Plasma) interactions, 
present in mid rapidity heavy ion collisions. 

It is important to note that at very forward rapidities {y > 3), one is probing the small 
X region of the nucleus and the large x region of the deuteron projectile. For example, at 
y ~ 3.2 and pt ~ 2, the deuteron wave function at x ~ 0.25 is probed. This is the region 
where valence quarks dominate over gluons and sea quarks in the deuteron wave function 
while the relevant x for the target nucleus is x ~ 4 x 10~^ so that gluons are the dominant 
parton species in the target nucleus. 

In this work, we concentrate in the very forward rapidity region {y > 3.2) so that 
the target nucleus is treated as a Color Glass Condensate while only the valence quarks 
in the projectile deuteron are included. We use the results of [7j for the scattering of 
valence quarks on a Color Glass Condensate to calculate charged hadron and neutral pion 
Pt spectra as well as the nuclear modification factor R^a in the forward rapidity region at 
RHIC. 

2 Scattering of quarks on a Color Glass Condensate 

In [7j, scattering of quark on a target described as a Color Glass Condensate was considered 
and the scattering cross section was calculated. The incoming quark is taken to be massless 
and carries zero transverse momentum. The scattering cross section is given by 
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Here, p~ is the light cone energy of the incoming quark while q" is the light cone energy 
of the outgoing quark with transverse momentum qt and adipoie is the cross section for 
scattering of a quark anti-quark dipole on a target described as a Color Glass Condensate. 
The dipole cross section satisfies the non-linear JIMWLK (BK |9J at large N^) evolution 
equation. To relate this quark-nucleus (proton) target scattering cross section to hadron 
production in deuteron gold collisions, we convolute this cross section with the quark 
distribution function in a deuteron and quark-hadron fragmentation function. We get 



where we have used the following kinematical relations Xq = k^e^ / z^/s, Xg = k^e~y / z^/s, 
Zmin = kx_ey 1^ (yi is the center of mass energy) while y and kt are the rapidity and 
transverse momentum of the measured hadron. Both the quark distribution and fragmen- 
tation functions depend on a factorization scale which is not written out explicitly. In 
our calculation, we set Q/ = kt where kt is the transverse momentum of the observed par- 
ticle. Eq. (j21) is the formula used in this work to calculate the hadron spectra at forward 
rapidities in deuteron gold collisions. 

We use the LO GRV98 quark distribution function JU] and the LO KKP quark hadron 
fragmentation function It should be noted that there is no available fragmentation 
function for negatively charged hadrons so therefore we use the fragmentation function for 
ih^ -|- h^)/2. However, in the transverse momentum range considered here {pt < 2GeV), 
this should not make a sizable difference. Since we are sensitive only to the large Xq 
(~ 0.2 — 0.5) region of the deuteron wave function, there is practically no sensitivity to the 
choice of the quark distribution function. Also, the nuclear modification of valence quarks 
in the deuteron wave function is minimal in this Xq range and is therefore neglected here. 

To proceed further, we need to know the dipole cross section adipoie which satisfies the 
JIMWLK equation. This is a very complicated functional equation ^2] which simplifies in 
the large A^^^ limit, known as the BK equation. The BK equation for the dipole cross section 
has been numerically solved by various people in different limits p!Bl. lancu, Itakura and 
Munier proposed a parameterization of the dipole cross section which has all the properties 
of the solution to the BK equation and used it to fit the HERA data on the proton structure 
function F2 ^14j. This is a very simple and economical parameterization (it basically has 
three free parameters, in addition to the light quark mass) which does an excellent job of 
describing the HERA data at x < 0.01. Therefore, we use this parameterization in this 
work. The dipole cross section is given by 
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The constants a, b are determined by matching the solutions at rtQs = 2 and js = 0.63 
and K = 9.9 are determined from LO BFKL. The form of the proton saturation scale Q^ 
is taken to be Ql = {xq/xY GeV^ with xq, A, A/q determined from fitting the HERA data 
on proton structure function F2. We refer the reader to for details of the fit. In case 
of a nucleus, we have Ra = 1.1 A^^^Rp and use the value of the minimum bias saturation 
scale given in [H] {Qlminbias — 0.95 GeV^^ at x = 0.01), with the same x dependence as in 
a proton given above. 

In Fig. (P), we show our results for charged hadron production aX y = 3.2 in proton 
proton collisions in arbitrary units while in Fig. Q, our results for charged hadron pro- 
duction in deuteron gold collisions is shown. Both figures are for minimum bias events. 
The proton-proton cross section is normalized to the data at fct = 1.1 GeV by multiplying 
hy & K factor of 2.57. To get the normalization of the dA spectra, we again multiply by 
a K factor which is slightly less than the = 1.8 factor for the proton-proton case. This 
is not unreasonable since Next to Leading Order corrections are typically large at low kt 
and that the K factors used can be different for nuclei ^Hl- However, the ratio RdA is 
calculated without any K factor. The agreement with the slope of the spectra for kf < 2 
GeV is quite reasonable specially since there are no free parameters in this calculation. 

In Fig. (jni, we show the nuclear modification factor RdA for hadron production in 
deuteron gold collisions at y = 3.2. We define 
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Again, the agreement with the data at low kt is quite good but the higher kt points start 
to show a deviation. This is discussed in more detail later. We emphasis the point that the 
ratio RdA is calculated without multiplying by any K factor and with no free parameters. 

It should be noted that the BRAHMS data is for negatively charged hadrons while 
our calculations are done for the average of positively and negatively charged hadrons 
since fragmentation functions for negatively charged hadrons are not available. Recently, 
Guzey et al. [12] investigated the dependence of the suppression and showed that isospin 
symmetry considerations can make a huge effect on the observed suppression. This effect 
is most pronounced in the kt > 2 GeV and can affect our results by (15 — 20)% in the 
kinematic region we cover. 

Finally, since the STAR collaboration has measured neutral pions at rapidity y = 3.8 
in deuteron gold collisions, we show our predictions for neutral pion nuclear modification 
factor RdA at ?/ = 3.8 in Fig. (jH). A slightly stronger suppression of RdA is seen at lowest 
kt as compared to charged hadrons at y = 3.2. 
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Figure 1: The invariant yield of charged hadrons at y = 3.2 at RHIC in proton-proton 
coUisions at ^/s = 200 GeV. The normahzation is a fit to the data. 
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Figure 2: The invariant yield of charged hadrons at y = 3.2 at RHIC in minimum bias 
deuteron-nucleus coUisions at y/s = 200 GeV. 



6 



en 



1.0 



0.8- 



< 

^^0.2 



0.0 




BRAHMS (h ) 



+ 



— CGC (h +h)/2 



1.0 



1.2 1.4 1.6 

(GeV) 



1.8 



Figure 3: The nuclear modification factor for cliarged hadrons at y = 3.2 in minimum bias 
deuteron-nucleus collisions at RHIC, y/s = 200 GeV. 
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Figure 4: The nuclear modification factor for neutral pions at y 
deuteron-nucleus collisions at RHIC, ^/s = 200 GeV. 
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3 Discussion 



In deuteron gold collisions in the forward rapidity and low kt region at RHIC considered 
in this work, the valence quarks are the most abundant parton species in the deuteron. 
They scatter on the target nucleus, which has its wave function fully developed (evolved 
in Xg as much as allowed by the kinematics) and is characterized by the nucleus saturation 
scale Qf{xg). The valence quark gets a transverse momentum kick of order Qf and is then 
"produced" . 

To get a rough idea of the scales involved and estimate where our approach should 
break down, we note that hadrons carry z fraction of the parent parton energy and that 
< z >~ 0.7 — 0.8 in this rapidity (for hadron production in pp collisions). The saturation 
scale of the nucleus is around Qf ~ 1.5 — 1.8 GeV (minimum bias) at rapidity of y = 3.2. 
The geometric scaling region [T71 HH] extends to a little higher momentum Qes = Ql/Qso 
and can be as high as 2.5 GeV. This means that our formalism should describe hadron 
production in minimum bias deuteron gold collisions in the forward rapidity region up to 
kt ^ z Qes ~ 2 GeV. We emphasize that these estimates are for minimum bias events only 
and the Color Glass Condensate formalism is expected to be valid at higher kt for more 
central collisions. 

As one goes to higher kt, gluon radiation becomes important jinj and will eventually 
dominate the hadron production cross section. This has not been included here and would 
presumably improve the high kt behavior of the spectra. Another caveat of our approach 
is the use of the dipole parameterization advocated in ^3]. This parameterization does 
not have the Cronin effect [2011^ 1221 • This may be partly responsible for the deviation of 
calculated RdA from the data. Also, this parameterization does not have the right high kt 
behavior since the double log limit is not built into it. This is the reason for the decrease of 
the RdA at higher kt which is not expected from pQCD. However, none of the other simple 
parameterizations of the dipole model j2S] which have been used to fit the HERA data have 
the right (BFKL) anomalous dimension. Therefore, we use this parameterization because 
it has the right anomalous dimension and we are staying in a limited kinematic region and 
since it is experimentally known that the Cronin effect goes away in the forward rapidity 
region. 

In order to extend this formalism to mid rapidity, one needs to include gluon production 
as considered in ^^1- The contribution of gluons to hadron production in mid rapidity is 
more important than the contribution of valence quarks since at mid rapidity, one probes 
small values of x in the deuteron wave function where there are a lot more gluons than 
quarks. Nevertheless, a practical problem with inclusion of gluons is that there is no 
available parameterization of the gluon-gluon dipole cross section which has been tested 
in other processes unlike the quark anti-quark dipole cross section which is probed in DIS 
processes. 

An important observable which has not been considered here is the centrality depen- 
dence of the hadron suppression factor in the forward rapidity region, recently shown by 
the BRAHMS collaboration in the Quark Matter 2004 pp. It will be very interesting to 
see whether our formalism can also describe the centrality dependence of the data. This 
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would require the application of the Monte Carlo Glauber method to our formalism since 
this is the method used by experimentalists to extract centrality dependence of the data. 
This is beyond the scope of this work and will be pursued later. 

The agreement of our calculations with the forward rapidity data becomes even more 
significant due to the fact we have not used any free parameters (for RaA and the slope 
of the spectra) in this calculation. All the necessary ingredients, the dipole cross section 
(for a proton) and the value of the nucleus saturation scale (for minimum bias events) 
have already been known for a while and used by various authors [HI E] ■ The fact that 
our simple and parameter free calculation based on the Color Glass Condensate formalism 
can reproduce the experimental data at low kt is a strong indication that the physics of 
forward rapidity region at RHIC is that of high gluon density QCD and the Color Glass 
Condensate. 
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